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We report structural and magnetic properties studies of large high quality single-crystals of the 
frustrated magnet, Nd 2 Zr 207 . Powder x-ray diffraction analysis confirms that Nd 2 Zr 207 adopts 
the pyrochlore structure. Room-temperature x-ray diffraction and time-of-flight neutron scattering 
experiments show that the crystals are stoichiometric in composition with no measurable site dis¬ 
order. The temperature dependence of the magnetic susceptibility shows no magnetic ordering at 
temperatures down to 0.5 K. Fits to the magnetic susceptibility data using a Curie-Weiss law reveal 
a ferromagnetic coupling between the Nd moments. Magnetization versus field measurements show 
a local Ising anisotropy along the (111) axes of the Nd 3+ ions in the ground state. Specific heat 
versus temperature measurements in zero applied magnetic field indicate the presence of a thermal 
anomaly below T ~ 7 K, but no evidence of magnetic ordering is observed down to 0.5 K. The 
experimental temperature dependence of the single-crystal bulk dc susceptibility and isothermal 
magnetization are analyzed using crystal field theory and the crystal field parameters and exchange 
coupling constants determined. 

PACS numbers: 61.50.-f, 75.10.Dg, 75.40.Cx 

PACS numbers: Valid PACS appear here 

I. INTRODUCTION 

Pyrochlore oxides of the type A% + B% + O-j (where 
A = trivalent rare earth, B = tetravalent transition metal 
element) are geometrically frustrated magnets^'-. These 
compounds have a face-centered cubic structure, with 
the space group Fd3m (No. 227). Both the rare earth 
atoms occupying the A sites and the transition metal el¬ 
ements located on the B sites form a pyrochlore lattice, 
a three-dimensional arrangement of corner-sharing tetra- 
hedra, known to display the highest degree of geometri¬ 
cal frustration, resulting in the many unusual magnetic 
properties of these systems. Depending on the nature of 
the interaction between the magnetic ions (A and/or B 
cations), these systems can exhibit either spin ice^, spin 
glass^, or highly correlated quantum disordered spin liq¬ 
uid states^, as well as long-range magnetic order—. 

One exciting avenue of research in this field has focused 
on systems where, due to a smaller magnetic moment 
on the trivalent rare earth A-site, quantum fluctuations 
play an important role in controlling the low temperature 
physics of the materials. Recent studies have pointed to 
the possibility of quantum spin liquid behavior (in which 


spin ice correlations exist at finite temperature, together 
with strong quantum effects), in a number of pyrochlore 
materials including Tb 2 Ti 20 7 and Yb 2 Ti 2 07 ^— . 

The success of research on pyrochlores and particu¬ 
larly the titanates, is due, in part, to the availability 
of large, high-quality single-crystals of these materials. 
Single-crystals of all the titanate pyrochlore family were 
successfully grown using the floating-zone technique^— 
and their magnetic ground states and properties have 
been investigated in great detail. (See Refs. [l| and [1 and 
references therein.) 

As the search for geometrically frustrated magnetic py¬ 
rochlores that exhibit quantum effects widens, the re¬ 
search community has turned its attention to the rare- 
earth zirconates, It has been shown that 

the rare-earth zirconates can be stabilized into two crys¬ 
tallographic structures, an ordered pyrochlore phase, for 
the first few elements of the lanthanide series (from lan¬ 
thanum to gadolinium), or a disordered fluorite phase, 
for the other lanthanide s 15 ’ 16 . In addition, these oxides 
undergo an order-disorder phase transition, from the py¬ 
rochlore to the fluorite structur e 15 ’ 17 . The temperature 
at which the structural transition occurs is strongly de- 
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pendent on the nature of the rare-earth element. In addi¬ 
tion, the co-existence of both the fluorite and pyrochlore 
phases has been observed in some zirconates with com¬ 
plex composition s 18 ' 19 . 

To date, due to the high melting point of the zir¬ 
conates, a majority of the studies of the structural and 
magnetic properties of the zirconate pyrochlores have 
been performed on powder samples— Recently, large 
single-crystals of the Pr containing zirconate pyrochlore, 
Pr 2 Zr 207 , have been grown by the floating-zone tech¬ 
nique^—. Studies of the magnetic properties have 
shown that the Pr 3+ ions in Pr 2 Zr 2 07 have a ground 
state doublet with a local Ising anisotropy along (111) 
axes and are coupled by an antiferromagnetic exchange, 
but no long-range magnetic order is observed down to 
very low temperatures (76 mK) 21 i 24 i 25 . The spin ice-likc 
correlations revealed by elastic magnetic neutron scat¬ 
tering, the broad excitation spectrum observed in the 
low-energy inelastic neutron scattering, and the reduc¬ 
tion of the nuclear contribution to the heat capacity^ 
may be connected with the strongly anisotropic superex¬ 
change and multipolar interaction s 26 ' 27 and interactions 
of the Pr 3+ ions with random lattice strains (strong sen¬ 
sitivity of the Pr 3+ ground state to lattice disorder in 
Pr 2 _ a; Bi a; Ru 207 pyrochlores was noted in Ref.Hsj). Some 
possible low temperature phases of the ProZmC^ py¬ 
rochlore have been discussed recently in Ref. 10. 

Neodymium zirconate, Nd 2 Zr 2 C> 7 , has been the sub¬ 
ject of just a small number of studies (see for example, 
Ref. Hd^. A study of the phase diagram of Nd 203 -ZrC >2 
shows that the pyrochlore oxide, Nd 2 Zr 2 07 melts congru- 
ently above 2000 °C, although the melting point was not 
established— The order to disorder transition in Nd 2 C> 3 - 
ZrC >2 occurs at 2300 ° C 15 ’ 3Q , raising the possibility that 
it may be difficult to prepare Nd 2 Zr 2 07 single-crystals 
of the pyrochlore structure due to the relatively small 
difference between the melting point and the structural 
transition temperature. When preparing this material it 
is therefore essential to establish the crystal structure of 
the Nd 2 Zr 2 0 7 samples and to adjust the conditions of 
synthesis in order to prepare the desired phase, as has 
been emphasized in previous report s 30 ' 31 . 

We have succeeded in preparing single-crystals of 
the neodymium zirconate pyrochlore, Nd 2 Zr 2 C >7 by the 
floating-zone technique. The availability of large high 
quality single-crystals is especially important for the 
study of the properties of this class of geometrically frus¬ 
trated magnets and, particularly, for solving the nature 
of the magnetic ground state. In this paper, we report 
the structural characterization and investigations of the 
magnetic properties of single-crystals of the neodymium 
zirconate pyrochlore, Nd 2 Zr 2 C> 7 . 


II. EXPERIMENTAL DETAILS 

Single-crystals of the Nd 2 Zr 2 07 pyrochlore oxide were 
prepared by the floating-zone technique, using a four- 


mirror Xenon arc lamp optical image furnace (CSI FZ-T- 
12000-X_VI-VP, Crystal Systems Incorporated, Japan). 
The growths were performed in air at ambient pressure 
and at growth speeds in the range 10-15 mrn/h. The 
crystal growth of Nd 2 Zr 2 07 is described elsewhere in a 
more detailed paper—. 

Room-temperature powder x-ray diffraction (XRD) 
measurements were performed on the as-grown boules in 
order to check the phase purity. Small quantities of the 
crystals were ground into powder and the data were then 
collected using a Panalytical x-ray diffractometer with 
Cu Kaq radiation (A = 1.5406 A), scanning 2d between 
10 and 110°, with a step size of 0.013° in 26 and a total 
scanning time of 16 hours. 

In order to investigate in detail the crystal structure, 
single-crystal diffraction experiments were carried out on 
small pieces cut from the as-grown boule of Nd 2 Zr 2 07 . 
Single-crystal x-ray diffraction measurements were per¬ 
formed at room temperature on an Oxford Diffraction 
Gemini diffractometer. A small piece of dimensions (0.28 
x 0.18 x 0.13 mm 3 ) was cleaved from the crystal boule 
and affixed to a glass fiber using epoxy resin. Data 
over large regions of reciprocal space were collected us¬ 
ing Mo Ka radiation (A = 0.71073 A). The single-crystal 
x-ray diffraction data were indexed and integrated us¬ 
ing CrysAlisPro (Agilent Technologies). The crystallo¬ 
graphic structure was refined using ShelXL—, as imple¬ 
mented in 01ex2— 

Neutron scattering measurements were performed us¬ 
ing the SXD beam line^, at the ISIS pulsed neutron 
source, Rutherford Appleton Laboratory. A crystal of 
cylindrical shape, 2.5 mm in radius and 5 mm in height, 
was cut from a Nd 2 Zr 2 C >7 crystal boule and mounted on 
an A1 pin with adhesive A1 tape and a piece of vanadium 
wire. The sample was placed in a closed-cycle helium 
refrigerator which was then evacuated to a vacuum of 
~ 10~ 5 mbar. Time-of-ffight neutron data were recorded 
at room temperature for a total of five orientations, with 
an exposure time of ~ 15 minutes per orientation. Data 
reduction was accomplished using the locally available 
software, SXD2001-. The structure was then refined 
using the ShelX software— 

Chemical composition analysis was carried out by en¬ 
ergy dispersive x-ray spectroscopy (EDAX) using a scan¬ 
ning electron microscope on two pieces cut from two ends 
of the Nd 2 Zr 2 07 crystal boule. 

The quality of the as-grown crystals was checked and 
the samples were aligned for cutting using a Laue x-ray 
imaging system with a Photonic-Science camera system. 
Rectangular prism shaped samples with the [100] (tetrag¬ 
onal), [110] (rhombic), and [111] (trigonal) directions per¬ 
pendicular to one face were cut from a Nd 2 Zr 2 C >7 single¬ 
crystal boule for the magnetization measurements. The 
demagnetizing factors for these rectangular prisms were 
calculated using expressions derived by Aharoni 36 . 

Magnetic susceptibility measurements as a function 
of temperature were carried out down to 0.5 K in ap¬ 
plied magnetic fields of 1 kOe using a Quantum De- 
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sign Magnetic Property Measurement System MPMS-5S 
SQUID magnetometer along with an i-Quantum 3 He in¬ 
sert. Magnetization measurements were also performed 
as a function of magnetic field up to 70 kOe directed 
along specific crystallographic axes at various tempera¬ 
tures. Heat capacity measurements in zero applied mag¬ 
netic field at temperatures from 0.5 to 400 K were carried 
out in a Quantum Design Physical Property Measure¬ 
ment System with a heat capacity option using a two-tau 
relaxation method. 


III. RESULTS AND DISCUSSION 
A. Single crystal growth and sample composition 

The crystal boules of the neodymium zirconate, 
Nd 2 Zr 2 0 7 , were transparent to light, with a dark pur¬ 
ple color. X-ray Laue diffraction patterns indicate that 
in most cases the facets of the crystals were aligned along 
the [Iff] direction. The large size and the good quality of 
the crystals make them suitable for most physical prop¬ 
erties characterization measurements, including those us¬ 
ing neutrons as a probe. 

Composition analysis by EDAX of the crystals shows 
that the cationic ratio averages close to 1:1 for Nd:Zr 
over the entire length of each crystal boule. The average 
atomic percentages of Nd, Zr, and O were (16.4 ± 0.3)%, 
(17.5 ±0.4)%, and (66.04 ±0.1)% respectively. Given the 
limitations of this technique there is reasonable agree¬ 
ment with the expected stoichiometry. 

B. Crystal structure 

Room-temperature x-ray diffraction data were col¬ 
lected on small pieces of crystals that were powdered. 
The x-ray data were fitted to the cubic FdSm space 
group^i, using the Rietveld refinement method, with the 
Fullprof software suited. The x-ray diffraction profile for 
the Nd 2 Zr 2 0 7 crystal displayed in Fig. |T] shows no trace 
of any impurity and contains a number of weak superlat¬ 
tice reflections demonstrating that the Fd3m pyrochlore 
structure is formed. These superlattice peaks, amongst 
which the strongest in the pattern are hkl = (111), 
(311), (331), and (511), indicate the cation and/or anion- 
vacancy ordering associated with the pyrochlore struc¬ 
ture^. Attempts to fit the x-ray data using different 
models for the occupancies of the Nd and the Zr sites, 
showed that the best model for the refinements is one in 
which there is no mixed occupancy of the cationic sites 
and no cationic/anionic deficiencies could be detected. 
Anisotropic atomic displacement parameters were em¬ 
ployed in the refinement of the crystal structure. The re¬ 
sults of the Rietveld refinement are given in Table [I] The 
lattice parameter was found to be slightly smaller than 
the previously reported values for polycrystalline sam¬ 
ples^. Analysis of the room-temperature x-ray diffrac- 
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FIG. 1. (Color online) Room-temperature powder x-ray 
diffraction pattern collected on a ground boule of Nd 2 Zr 2 0 7 . 
The experimental profile (red circles) and a Rietveld refine¬ 
ment (black line) made using the cubic Fd3m structure are 
shown, with the difference given by the blue line. The posi¬ 
tions of the Bragg peaks are indicated by green vertical bars. 
The positions of the (111), (311), (331), and (511) pyrochlore 
superlattice reflections are indicated by the violet arrows. 


TABLE I. Crystallographic parameters for Nd 2 Zr 2 0 7 refined 
using room-temperature powder x-ray diffraction data. 


Structural data 

Crystal system 

Cubic 

Space group 

Fd3m 

a (A) 

10.62657(7) 

Agreement factors 

R P 

7.99% 

R-wp 

10.5% 

Rex. p 

6.13% 

s 

1.71 


tion data collected on the powder used to prepare the 
feed rods showed that the calculated value of the lattice 
parameter (a = 10.68287(5) A) is, in this case, very sim¬ 
ilar to that reported in the literatur e 29 ! 37 . 

The powder x-ray diffraction data measurements are 
sensitive to any preferred orientation present in the sam¬ 
ple, and this is reflected in changes in the relative inten¬ 
sities of some of the Bragg peaks. In order to confirm the 
crystallographic parameters obtained from the powder x- 
ray diffraction profile, a small piece was cleaved from the 
crystal boule and single-crystal x-ray diffraction measure¬ 
ments were performed at room temperature. Attempts 
to fit the single-crystal x-ray data using different mod¬ 
els for the occupancies of the cations sites showed that 
there is no deficiency in Nd or Zr content, and no internal 
site disorder was detected in our Nd 2 Zr 2 0 7 crystals. The 
best model for the refinement of the occupancies of the 
anion sites was found to be one in which only the O 48/ 
site occupancy varies. Anisotropic thermal displacement 
parameters were used to refine the crystal structure. The 
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TABLE II. Refined structural parameters for Nd 2 Zr 2 C >7 
from room-temperature single-crystal x-ray diffraction data. 
Atomic positions for the Fd3m (origin choice 2) cubic struc¬ 
ture are Nd, 16d (|, |); Zr, 16c (0, 0, 0); O, 48/ (x, |, |); 

and O', 8b (|, |, |). Atomic displacement parameters in units 

(1(T 3 A 2 ). 


Atom 

U iso / U eq 

Uu 

U22 = U33 U12 = U13 

U23 

Nd 

5.61(15) 

5.61(15) 

5.61(15) 

-1.72(3) 

-1.72(3) 

Zr 

6.4(2) 

6.4(2) 

6.4(2) 

2.50(6) 

2.50(6) 

O 

11.6(4) 

15.5(10) 

9.7(5) 

0 

2.4(8) 

O' 

8.5(7) 

8.5(7) 

8.5(7) 

0 

0 


a = 10.62652(9) A, R = 2.32%, = 1-224. 

O SOF = 0.936(8), Oi = 0.3360(2). 


TABLE III. Refined structural parameters for Nd 2 Zi '207 from 
room-temperature neutron diffraction datcAA Atomic posi¬ 
tions for the Fd3m (origin choice 2) are the same as in Ta¬ 
bic EH Atomic displacement parameters in units (10 3 A 2 ). 


Atom 

Uiso / Ueq 

U u 

U22 — U33 U\2 — u 13 

U23 

Nd 

5.6(5) 

5.6(5) 

5.6(5) 

-2.0(5) 

-2.0(5) 

Zr 

9.3(6) 

9.3(6) 

9.3(6) 

2.8(7) 

2.8(7) 

O 

11.5(6) 

15.5(12) 

9.6(7) 

0 

3.1(9) 

O' 

6.9(9) 

6.9(9) 

6.9(9) 

0 

0 

a = 10. 

.618(3) A, 

R = 6.93%, x* = : 

1.627. 


O SOF 

1 = 0.97(1), Oi = 

0.3360(3) 




structural parameters obtained from the best fit of the 
x-ray diffraction data at room temperature are listed in 
Table EU 

In order to confirm the structural model obtained using 
the single-crystal x-ray diffraction data, a small piece of 
an as-grown single-crystal of Nd 2 Zr 2 C >7 was cut from the 
boule and single-crystal time-of-flight neutron diffraction 
data were collected at room temperature. Attempts to fit 
the neutron data using different models for the occupan¬ 
cies of the Nd and Zr sites, showed that the best model 
for the refinements is one in which the cationic occupan¬ 
cies are fixed to their nominal value of 1. (When the site 
occupancy factors (SOF) were allowed to vary, the re¬ 
finement revealed an over-occupancy of the two cationic 
sites.) The occupancy of the oxygen 8b site was fixed to 
the value of 1, as the oxygen vacancies are usually found 
to be mainly on the oxygen 48/ site^S. (This was pre¬ 
viously confirmed by the refinement of the single-crystal 
x-ray diffraction data.) The refined value of the site oc¬ 
cupancy factor for the 48/ site occupied by the oxygen 
shows that there are very few oxygen vacancies in the 
crystallographic structure of our Nd 2 Zr 2 C >7 crystals. The 
crystal structure was refined using anisotropic thermal 
displacement parameters. A summary of the results of 
the Rietveld refinement of the neutron diffraction data 
is given in Table Mil The structural parameters obtained 
using the neutron diffraction data were found to be con¬ 
sistent with those calculated using x-ray diffraction data. 

The value of x = 0.3360(3) for the O atom in the 48/ 
site confirms that both the polyhedra formed around the 


A (Nd) and B (Zr) sites are distorted from the ideal ge¬ 
ometries (a perfect cube and a perfect octahedron respec¬ 
tively)^. This fact is strongly reflected in the uniaxial 
nature of the A site symmetry and in the magnetic prop¬ 
erties of the system. The scalenohedron NdOs is formed 
by two short bonds Nd-O' of 2.2988(2) A length, and six 
bonds Nd-0 of 2.561(2) A, while the trigonal antiprism 
ZrC>6 has all six bonds Zr-0 of 2.087(1) A. The two short 
bonds Nd-O' are stacked along the [111] direction, the 
other six long bonds (Nd-O) connect the Nd ion at the 
center with oxygen (O) ions at the apexes of the trigonal 
antiprism strongly compressed along the [111] axis. 

The structural features of the Nd 2 Zr 2 07 crystals are 
very similar to those of the neighboring praseodymium 
zirconate, Pi' 2 Zr 2 07 . Rietveld refinement of the neutron 
diffraction data collected at the SXD beam line on a 
Pr 2 Zr 2 07 single-crystal indicates a value of x = 0.3349(1) 
for the O atom in the 48/ siteAL. The oxygen polyhedra 
formed around the Pr and the Zr ions show a similar de¬ 
gree of distortion to the those of neodymium zirconate. 
The bond lengths inside the two coordination polyhe¬ 
dra are as follows: 2.3088(2) A for the Pr-O' bonds, 
2.578(1) A for Pr-0 bonds, and 2.0917(8) A for the Zr-0 
bonds. Given the similarities between the structural pa¬ 
rameters of the Nd 2 Zr 207 and Pr 2 Zr 207 pyrochlores, we 
can expect some similarities between the crystal fields in 
these two systems. 

In summary, Rietveld refinements of the powder and 
single-crystal x-ray diffraction and neutron diffraction 
data confirm that Nd 2 Zr 2 07 adopts the cubic Fd3m 
pyrochlore structure. The results show no measurable 
cationic or anionic deficiencies of the as-grown crystals 
of Nd 2 Zr 2 0 7 and confirm the good quality of the boules 
and their suitability for future investigations, including 
those using neutrons as a probe, given the large size of 
the crystal boules grown by the floating-zone technique. 


C. Magnetization and heat capacity 

Field-cooled (FC) and zero-field-cooled (ZFC) mag¬ 
netization versus temperature curves were measured on 
a Nd 2 Zr 2 C >7 single-crystal aligned along three different 
crystallographic directions ([100], [110], and [111])- The 
temperature dependence of the dc magnetic suscepti¬ 
bility, x(^)> an d reciprocal dc magnetic susceptibility 
X” 1 (T) are shown in Fig. [2j No anomalies in x(^) are 
observed down to 0.5 K suggesting the absence of a mag¬ 
netic transition. The magnetic susceptibility data mea¬ 
sured along the different directions all exhibit a mono¬ 
tonic decrease when warming from 0.5 to 350 K, and at 
temperatures above 5 K the magnetic susceptibilities col¬ 
lected in 1 kOe along the three high symmetry directions 
all overlap to within experimental error. 

Attempts to fit the temperature dependent reciprocal 
magnetic susceptibilities reveal that the x -1 (T) data do 
not obey a Curie-Weiss law in the temperature range 
0.5 to 350 K. Nevertheless, we have made fits to a Curie- 
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FIG. 2. (Color online) (a) Temperature dependence of the dc magnetic susceptibility, x versus T, in the temperature range 
1.8 to 350 K for a crystal of Nd 2 Zr 2 07 , with a magnetic field applied along the [100] (red), [110] (black), and [111] (orange) 
directions. The inset shows x versus T in the temperature range 0.5 to 1.75 K. (b) Measured (symbols) and calculated (the solid 
curve) temperature dependence of the reciprocal of the bulk dc susceptibility, x _1 versus T, of Nd 2 Zr 207 for a field applied 
along the [110] direction. The temperature dependence of the reciprocal bulk dc susceptibility was calculated using the crystal 
field parameters determined in the present work. The inset shows x _1 versus T and the linear fit (using the Curie-Weiss law) 
to the data in the temperature range 1.8 to 10 K for a magnetic field applied along the [111] direction. 


Weiss law over a reduced temperature range (1.8 to 10 K) 
(see Fig.d> inset), and the results of these fits have shown 
that for T < 10 K Nd 2 Zr 2 07 has an effective moment of 
/r e ff = 2.543(2 )/tb (^b is the Bohr magneton) and a Weiss 
temperature of f?w = +0.200(8) K. Over an extended 
temperature range from 1.8 to 60 K, the measured tem¬ 
perature dependence of the dc susceptibility can be ap¬ 
proximated by the expression \ ( T ) = C/ (T — $w)+Xvv 
with the Weiss temperature of $w = +0.150 K and 
a Van Vleck contribution of \vv = 0.00346 emu/(mol 
Nd). The effective moment of Nd 3+ in Nd 2 Zr 2 07 is 
smaller than the magnetic moment of -^=/iB of a free 

Nd 3+ in the ground state 4 /g/ 2 . The $w value indi¬ 
cates a ferromagnetic coupling between the Nd spins. 
We have also performed magnetic susceptibility measure¬ 
ments on a ground piece of the Nd 2 Zr 2 07 crystal boule. 
The estimated Weiss temperature 0\y was, in this case, 
found to be +0.108(6) K. Previous studies also reported 
small absolute values for the Weiss temperature in Nd- 
based pyrochlores, with dyv = +0.06 K in Nd 2 Zr 2 072 £, 
6> w = -0.069(4) K in Nd 2 Pb 2 0 A and 6> w = -0.31 K in 
Nd 2 Sn 2 0 7 43 i 44 . Small differences in the absolute values 
obtained for Nd 2 Zr 2 07 may be explained by the fact that 
some measurements were carried out on aligned crys¬ 
tal, along well-defined crystallographic axes, while others 
were for polycrystalline samples. The results also appear 
to depend on the exact temperature range over which the 
fits are performed. The difference between the data col¬ 
lected for different directions of the applied field at low 
temperatures (see inset in Fig. H-) cannot be explained by 
the contribution of the demagnetizing field; the difference 
in the susceptibilities corrected for the demagnetization 


factor has been found to be close to 1%, for all the crys¬ 
tallographic directions. However, this difference, which 
is found in a small but finite magnetic field of 1 kOe, 
and which increases with decreasing temperature, may 
signal the onset of a phase transition at lower tempera¬ 
ture. Different signs of $w hr compounds with different 
S-site cations (a chemical pressure effect) may be caused 
by closely competing exchange and dipolar interactions 
between the Nd 3+ ions. 

Figure [3] shows the magnetization as a function of ap¬ 
plied magnetic field M ( H ) at a temperature of T = 0.5 K 
for Nd 2 Zr 2 07 along three high symmetry directions, 
[100], [110], and [111]. The magnetization response is 
reversible with no hysteresis between the field increasing 
and field decreasing M (H) curves. The data reveal a 
non-linear variation of the magnetization as a function 
of applied Held. The magnetization is highest in strong 
magnetic fields applied along the [100] direction and low¬ 
est for the [110] direction. The field dependence of the 
magnetization along the three high symmetry directions, 
measured at various temperatures (see Fig. 11 indicate 
that the magnetic anisotropy evolves on cooling. The 
results of the magnetization measurements suggest that 
the magnetization follows a local Ising behavior, with a 
different dependence of M ( H ) along different crystallo¬ 
graphic directions. This kind of magnetic anisotropy rep¬ 
resents a key feature of the spin ice system, as observed 
in the Dy 2 Ti 2 07 and Ho 2 Ti 2 07 pyrochlore oxide s 45 ’ 46 . 
However, the field dependent magnetization data reveal 
that the magnetic anisotropy observed for the Nd 2 Zr 2 07 
pyrochlore is only partially consistent with the spin ice 
model, as the values of the magnetic moments measured 
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FIG. 3. (Color online) Isothermal magnetization (M) as a 
function of applied magnetic field (H) along the [100] (red), 
[110] (black), and [111] (orange) directions at T = 0.5 K for 
a single-crystal of Nd 2 Zr 2 07 . 


for the [100], [110], and [111] directions, at the maximum 
applied field do not agree exactly with the ratios of the 
expected saturated moments for a classic spin ice configu¬ 
ration^ (see Fig. 0K) and that in contrast to the titanate 
spin ice materials, there is no evidence for a shoulder de¬ 
veloping in M (H) when H is applied along [111]. These 
results are, however, similar to those reported for the re¬ 
lated system Pr 2 Zr 2 07 22 ’ 24,2 - and provide evidence for 
a mixing of the wave functions of the ground-state dou¬ 
blet and excited crystal field (CF) levels of the Nd 3+ and 

p r 3 + 

ions in strong external magnetic fields. 

Specific heat measurements were performed in the tem¬ 
perature range 0.5 to 400 K. The temperature depen¬ 
dence of the heat capacity of Nd 2 Zr 2 07 in zero external 
magnetic field is shown in Fig. [5j The data agree well 
with the previously published results for Nd 2 Zr 2 O y 20 ' 47 . 
A continuous decrease in the specific heat with decreas¬ 
ing temperature is followed by an upturn at temperatures 
below ~ 7 K. Previous heat capacity studies performed 
on a polycrystalline sample of Nd 2 Zr 2 C >7 show a phase 
transition centered at 0.37 K and evidence of some mag¬ 
netic ordering— 

D. Crystal field parameters 

To analyze the magnetic behavior of Nd 2 Zr 2 07 , we 
consider the following Hamiltonian of a single Nd 3+ ion 
in an applied magnetic field H: 

h = tipi + hep + hz- (1) 

Here hpi is a free ion standard Hamiltonian^ that oper¬ 



tf(kOe) //(kOe) 


FIG. 4. (Color online) Measured (symbols) and calculated 
(solid lines) external magnetic field dependence of the magne¬ 
tization (M) of Nd 2 Zr 2 0 7 at various temperatures. The field 
was applied along the [100] (red), [110] (black), and [111] (or¬ 
ange) axes of the Nd 2 Zr 2 07 crystal lattice, at temperatures 
of (a) 1.8 K, (b) 5 K, (c) 10 K , and (d) 20 K. The field depen¬ 
dence of the magnetization was calculated using the crystal 
field parameters determined in the present work. The dashed 
lines in (a) indicate the theoretically expected values of the 
saturated magnetization satisfying the ratios v / 2:^/3/2:l for 
the [100], [111], and [110] directions, respectively, in a classical 
spin ice configuration— 


ates in the total space of 364 states of the electronic 4/ 3 
configuration. The next term in ©, 

hcF = B° 2 C 2 0 + H 4 °C 4 + B\(Ci - Cl 3 ) + H 6 °C 0 6 

+bUc!-c 6 _ 3 ) + bI(c 6 6 + c 6 _ 6 ), (2) 

where hcF is the energy of 4/ electrons in the crystal 
field of the perfect lattice with the D 3 d symmetry at the 
Nd 3+ sites written in the local system of coordinates with 
the z axis along the trigonal symmetry axis. B£ are CF 
parameters and C £ are spherical tensor operators. The 
third term in © is, 

hz = /z ■ ffioc. (3) 

where hz is the electronic Zeeman energy. fi = 
—/xb (kL + 2 S) is the magnetic moment of a Nd 3+ ion 
(L and S are the electronic orbital and spin moments, 
respectively, and k is the orbital reduction factor—). The 
local magnetic held affecting the Nd 3+ ions, H\ oc = 
H + H ex ch + Hdip — -Hp, involves the external held, H , 
the exchange and dipolar helds, H ex c h and iTdipi respec¬ 
tively, corresponding to anisotropic exchange and mag¬ 
netic dipolar interactions between the Nd 3+ ions that are 
considered within the self-consistent held approximation 
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FIG. 5. (Color online) Heat capacity as a function of tem¬ 
perature in zero applied magnetic field for a Nd 2 Zr 20 r single 
crystal. The inset shows the temperature dependence of the 
heat capacity data on a logarithmic scale. 


, and the demagnetizing field H d = AnNM/3 , where N 
is the demagnetizing factor, M = Y2 n m n/v is the mag¬ 
netization, and v = a 3 /4 is the unit cell volume. There 
are four magnetically non-equivalent rare-earth ions in 
the unit cell and m n is the average magnetic moment of 
a Nd 3+ ion belonging to the sublattice n. 

The dipolar fields at sites n, Hdip.na = 
J2n>p PQn£ m n'0 (P = Air/3v), are determined by 

dimensionless lattice sums which have been com¬ 

puted in Ref. |4^. The exchange interaction is assumed 
to be non-zero for the nearest neighbor Nd 3+ ions 
only. In particular, the exchange interaction between 
the ions with radius vectors rq = (|,A,A)a and 

r 2 = (—g, — |, |) a is approximated by 

^exch(f;2) = \\Plx'Plx' ^-LlPlz 1 P2z' 

-A_L 2 Ml 3 /'M 2 y>, (4) 

where p x >, p y ', and p z ' are operators of the magnetic 
moment components along the vectors X' = rq — r 2 , 
Y' = r 2 xri, and Z' = r i + rq. To reduce the number 
of unknown parameters, we neglect in ([4]) an additional 
antisymmetric term Adm^ 7 • (Mi x M 2 ) (Dzyaloshinskii- 
Moriya like) which is allowed by symmetry of the py- 
rochlore lattice^. The values of the exchange coupling 
constants Ay, A_li, and Aj_ 2 , six CF parameters, and 
the orbital reduction factor were obtained from the fit¬ 
ting procedure by comparing the measured and com¬ 
puted magnetic field dependence of the heat capacity, the 
bulk dc susceptibility with temperature, and the mag¬ 
netic field dependence of the isothermal magnetization. 

Modeling of the magnetization, dc susceptibility, and 
the contribution of the Nd subsystem to the heat capac¬ 



FIG. 6. Calculated and experimentally estimated contribu¬ 
tions to the heat capacity of Nd 2 Zi’ 207 , due to CF excitations 
in the Nd 3+ ions ions (the temperature range is extended 
up to 650 K using the data taken from Ref. 154] '). The solid 
and dashed curves correspond to the results of calculations 
based on the CF energies obtained in the present work and 
in Ref. |56|, respectively. 


ity involved numerical diagonalization of the Hamilto¬ 
nian © for fixed values of the applied magnetic field and 
subsequent quantum-statistical averaging of the single¬ 
ion energy, <h>, the square of energy, <h 2 >, and 
the magnetic moment components for different temper¬ 
atures. Parameters of the free ion Hamiltonian, hp 1 , 
were taken from Ref. [5l]. The initial values of the CF 
parameters were calculated in the framework of the ex¬ 
change charge mode l 52 i 53 by making use of the parame¬ 
ters of the crystal lattice determined in the present work 
(see column A in Tablg JIAT . The effective ionic charges 
were taken from Ref. [53| and the scaling factors which 
determine the exchange charges on the Nd-O' and Nd- 
O bonds were obtained from the analysis of the con¬ 
tribution to the specific heat from the CF excitations 
in the Nd subsystem. This contribution as a function 
of temperature was estimated from the difference be¬ 
tween the measured heat capacities of Nd 2 Zr 2 07 and 
non-magnetic La2Zr20y given in Ref. HH- The results 
of calculations of the heat capacity of the Nd subsys¬ 
tem Osid = -Na ( <h 2 > — <h> 2 ) fkp,T 2 ( Na is the Avo- 
gadro number, is the Boltzmann constant) match well 
the experimental data for temperatures below 300 K (see 
Fig. 0. More detailed analysis, in particular taking into 
account changes of phonon frequencies^, is necessary to 
clarify the behavior of the relatively small differences be¬ 
tween the heat capacities of neodymium and lanthanum 
zirconates at higher temperatures. 

The final values of the CF parameters which were used 
in our simulations of the dc susceptibility and the magne¬ 
tization for different directions of the external magnetic 
fields are presented in Table EO column B. These pa- 








TABLE IV. Crystal field parameters Bp (cm -1 ) in some Nd 
and Pr based pyrochlores. The initial CF parameters for 
Nd 2 Zr 2 07 , calculated using the exchange charge mode l 52 ' 53 
with the parameters of the crystal lattice determined in the 
present work, are given in column A. The final values of the 
CF parameters used to simulate the dc susceptibility and the 
magnetization data are given in column B (see text for de¬ 
tails). 


p q Pr 2 Sn 207 Pr 2 Zr 207 Nd 2 lr 2 07 Nd2Zr2C>7 
Ref. 57 Ref. 24 Ref. 58 Present work 

A B 


2 0 

466.3 

711 

702 

438 

460 

4 0 

3485.5 

3173 

2110 

3096 

3200 

4 3 

1296.6 

1093 

3087 

1027 

1132.5 

6 0 

1163.7 

1402 

1500 

1053 

992 

6 3 

-865.7 

235 

-284 

-608 

-593.4 

6 6 

1548.7 

604 

1007 

804 

1021 


rameters agree qualitatively with those found in Ref. [5?] 
from the analysis of the CF excitations in the Pr 3+ ions 
in Pr 2 Sn 2 07 but differ markedly from the results of the 
studies of the crystal fields in Nd 2 Ir 2 0v^ (as far as we 
know, this is the only earlier published set of CF pa¬ 
rameters for Nd 3+ ions in Nd-based pyrochlores) and 
Pr 2 Zr 2 0v^ (see Table HVl) . The CF interaction has been 
treated in Refs. and i58l as a perturbation within the 
ground J-multiplet only. However, the CF splittings of 
the ground multiplets 4 /g/ 2 of Nd 3+ and 3 I7 4 of Pr 3+ 
ions in the pyrochlores are rather large and comparable 
to the excitation energies of the nearest multiplets 4 /n/ 2 
and 3 H 5 , respectively. As a result, the J-mixing should 
be taken into account. 

The calculated energies of the sub-levels of the 4 /g/ 2 
ground multiplet are 0 (r 56 ), 178.5 (r 4 ), 252 (r 56 ), 262 
(r 4 ) and 785 (r 4 ) (in cm -1 ). The corresponding irre¬ 
ducible representations of the D%d point symmetry group 
are presented in brackets. A transverse g-factor of the 
Kramers doublet with the Tsg symmetry (of dipolar- 
octupolar type according to the classification proposed 
in Ref. [59f ) is exactly zero and our results are in line with 
the Ising-type magnetic anisotropy of the Nd 3+ ions in 
compounds with the pyrochlore structure noted previ¬ 
ously in Ref. [6(| The calculated longitudinal ^-factor of 
the ground state, g\\ = 4.793, is close to the value 4.6 
that has been determined in Ref. [fil] for the Nd 3+ ions in 
Nd 2 Mo 2 07 . 

The temperature dependence of the dc susceptibility 
was calculated using the corresponding expression pre¬ 
sented in Ref. [49|, the results of the calculations are shown 
in Fig. [2jo. In order to match the experimental data, 
the orbital reduction factor k = 0.9575 and the exchange 
constants Ay = 0.2, A_li = -0.1 and Aj_ 2 = 0 (in units of 
kO e/ns) were introduced. 

The calculated Van Vleck susceptibility of the Nd 3+ 
ions in the ground state, averaged over four magneti¬ 
cally non-equivalent sites in the pyrochlore lattice, xvv = 
2 T,j, a =x, v ,z 10 ' I Me* | g)\ 2 /3 (.Ej - E g ) (here Ej and E g are 


the energies of the excited states (j) and the ground state 
(g)), xvv = 0.00327 emu/(mol Nd), matches very well 
with the value obtained from the analysis of the experi¬ 
mental data (see Fig. [2] and subsection IIII Cl above). 

Because a gap of 250 K between the ground and the 
first excited doublets in the energy spectrum of Nd 3+ 
ions is rather large, the low temperature magnetic prop¬ 
erties of the Nd subsystem in Nd 2 Zr 2 C >7 can be consid¬ 
ered using the truncated Hilbert space spanned by the 
wave functions of the ground doublets only. The projec¬ 
tion of the exchange and magnetic dipolar interactions 
between the nearest neighbor Nd 3+ ions on this space 
can be written in terms of the pseudospin S = 1/2 op¬ 
erators defined in local Cartesian coordinates with the 
z-axes along the corresponding trigonal symmetry axes 
crossing at the angle of arccos(—1/3): 

H s = (Jexch + Jdd) &zlS z 2, (5) 

heie J exch — ( 2 ^ll Aj_ 4 ) /3 and Jdd — 

l r i — T 2 \ ■ Using the values of parameters given 
above, we obtain J exc h = 0.257 K and J^d = 0.442 K. 
So, both the exchange and dipolar interactions be¬ 
tween the nearest neighbor ions have antiferromagnetic 
character, and the exchange contribution is about half 
the dipolar contribution. At the first glance this re¬ 
sult contradicts with the sign of the Weiss tempera¬ 
ture indicating a ferromagnetic coupling. However, in 
case of Ising-type magnetic anisotropy along the crys¬ 
tal ternary axes in the pyrochlore structure, the sum 
of antiferromagnetic interactions between a fixed rare- 
earth ion and the six nearest neighbor ions leads to for¬ 
mation of a local magnetic field parallel to the applied 
field and, correspondingly, to a positive Weiss temper¬ 
ature^. A crude estimate, taking into account only 
the ground doublet of the Nd 3+ ions and neglecting 
the Van Vleck contribution to the susceptibility, gives a 

value $w = (s^Mb/^b) [2A|| — A_li + 7.52 (47t/ 3^)] = 
0.38 K with the dominant contribution (0.25 K) from 
the dipolar interactions represented by the last term in 
the square brackets. According to a general form of the 
exchange interaction^, it may contain terms (neglected 
in the present work) constructed from single-ion opera- 

tors E/C'l (/)«+(/), E/ C -2 (/)«-(/)> E/Cf (/)*+(/), 

E/ U® 5 (/)s_ (/) (the summation over / is over individ¬ 
ual 4/ electrons with spin s (/) on the Nd 3+ ion) which 
mix eigenfunctions of the pseudospin S z operator— and 
bring about fluctuations of the Ising-type magnetic mo¬ 
ments. The effects of these terms on the low tempera¬ 
ture magnetic properties of the Nd 3+ ions remain to be 
explored. It should also be noted that at temperatures 
below 1.5 K the single-site self-consistent approximation 
appears not to be valid (in particular, in contrast to the 
experimental data which demonstrate \ (T) doubles in 
the temperature range from 1 to 0.5 K, the calculated 
value of x (T) increases more rapidly as the temperature 
decreases below 1 K). 









9 


The simulated magnetic field dependence of the mag¬ 
netization of Nd 2 Zr 2 07 with the magnetic field directed 
along the [100], [110], and [111] symmetry axes of the py- 
rochlore lattice at different temperatures are compared 
with the experimental data in Fig. [4] As can be seen 
from Fig. [4) the calculated field dependences of the mag¬ 
netization agree well with the experimental data. 


IV. SUMMARY 

We have investigated the structural properties of 
single-crystals of the pyrochlore oxide, Nd 2 Zr 2 07 . Pow¬ 
der x-ray diffraction studies of the crystal boule confirm 
that the Fd3m pyrochlore phase is formed, and single¬ 
crystal x-ray diffraction and neutron diffraction data 
show a Nd 2 Zr 2 07 structure with no cationic or anionic 
deficiencies. Supplementary investigations, such as de¬ 
tailed inelastic neutron scattering experiments, are now 
being carried out on these Nd 2 Zr 2 07 crystals. 

The magnetic susceptibility data collected reveal a spin 
ice-type magnetic anisotropy and an effective ferromag¬ 
netic coupling between the Nd spins. The heat capac¬ 
ity decreases monotonically with decreasing temperature, 
followed by an upturn at low temperature. However, no 


sign of long-range magnetic ordering was observed in ei¬ 
ther the heat capacity or the magnetization data down 
to 0.5 K. 

The measured temperature dependence of the mag¬ 
netic susceptibility and the field dependences of the 
isothermal magnetization are reproduced satisfactorily 
by the CF calculations and subsequent modeling in the 
framework of the mean-field approximation taking into 
account the dipolar and the anisotropic exchange inter¬ 
actions between the Nd 3+ ions. The sets of CF parame¬ 
ters and exchange coupling constants obtained may serve 
as a basis for future studies of the spectral and magnetic 
properties of Nd-based pyrochlore magnets at low tem¬ 
peratures. 
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